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ABSTRACT 

We present the results of an exploratory study of broad line region (BLR) metallicity in 34 
2.2 ^ z ^ 4.6 quasars with far-infrared (FIR) luminosities (Lfir) from lO^^ '* to ^lO'^'^ L©. 
Quasar samples sorted by LpjR might represent an evolutionary sequence if the star formation 
rates (SFRs) in quasar hosts generally diminish across quasar lifetimes. We use rest-frame 
ultraviolet spectra from the Sloan Digital Sky Survey to construct three composite spectra 
sorted by LpiR, corresponding to average SFRs of 4980, 2130 and ^ 340 M© yr^^ after cor- 
recting for a nominal quasar FIR contribution. The measured N V A1240/C IV A 1550 and 
Si IV A 1397H-0 IV] A 1402/C IV A 1550 emission Hne ratios indicate super-solar BLR metal- 
licities in all three composites, with no evidence for a trend with the star formation rate. The 
formal derived metallicities, 5-9 Z©, are similar to those derived for the BLRs of other 
quasars at similar redshifts and luminosities. These results suggest that the ongoing star for- 
mation in the host is not responsible for the metal enrichment of the BLR gas. Instead, the 
BLR gas must have been enriched before the visible quasar phase. These results for high 
quasar metallicities, regardless of LpiR, are consistent with evolution scenarios wherein visi- 
bly bright quasars appear after the main episode(s) of star formation and metal enrichment in 
the host galaxies. Finally, young quasars, those more closely associated with a recent merger 
or a blowout of gas and dust, may exhibit tracers of these events, such as redder continuum 
slopes and higher incidence of narrow absorption lines. With the caveat of small sample sizes, 
we find no relation between Lfir and the reddening or the incidence of absorption lines. 
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1 INTRODUCTION 

Numerous studies have shown a strong correlation between su- 
permassive black hole (SMBH) mass and host galaxy mass over 
broad mass ranges, implying that there is a close evolutio nary rela- 
tionship between them (Xjebhardt et al. 2000; Merritt & Ferraresd 



2001 



lTremaineetai]|2002l : iMarconi & Hunlll2003l : iHaring & Rij3 



2004 [Shields et al.l 1200^ . Luminous quasars at high redshifts. 



which represent the most massive SMBHs, are experiencing vigor- 
ous accretion of a significant portion of the final SMBH mass. This 
accretion activity is believed to be triggered by global processes 
(galaxy mergers, interactions or perhaps secular evolution) that also 
trigger major episodes of star formation in the massive host galax- 
ies, which are also rapidly being assembled at high redshifts. There- 
fore, the processes of SMBH formation and growth resulting in 
the quas ar phenomenon are directly linked to the birth of massive 
galaxies iHaehnelt et alJI 19981; iRichstone et al][ 19981 : lOmont et all 
l200lLl2003l ; lBeelen et alJl2006l ; ICox et al.ll2006l) . However, the na- 
ture of the relationship between quasars and galaxy formation is 
not well understood. 

It is widely believed that strong interactions between gas rich 



* email: 
(FH) 



lsiinon9astro.ufl.edu (LES); hamaiin@astro.ufl.edu 



galaxies can result in ultra-luminous infrared galaxies (ULIRGs), 
defined by Lir >10' ^L©, which have star formation rates ( SFRs) 



of > 100 Mr:, yr 
Flores etal..,2004; 



(Houcketal. 1985 
Cox et al. 2005; Beelen et al 



Omont et al.'"200l',%0: 



2006;.Daddiet 



I 



2007l ; ICaoetalJl2008l) . A fraction of ULIRGs have been found 



to contain dust-enshrouded quasars dSanders & Mirabell 1 19961 ; 
I Lonsdale^ al. 2006). There is evidence at both low and high red- 
shift that these embedded q uasars are precu rsors to optically lumi- 
nous quasars. For example. ICao et all ( I2OO81) find that low redshift 
quasars with IR luminosities in the ULIRG range fall between op- 
tically selected PG quasars and ULIRGs (starbursts) for a variety 
of mid-IR spectroscopic indicators of starburst and active galactic 
nucleus (AGN) contributions including: polycyclic aromatic hydro- 
carbon (PAH) luminosities, fine structure emission line strengths, 
silicate absorption strengths, spectral slope and mid-IR colour in- 
dices. These findings suggest that IR luminous quasars could rep- 
resent an intermediate stage between a dominant starburst and a 
dominant AGN phase. 

Quasars with high SFRs therefore may be at an intermediate 
stage between an embedded quasar in a star forming galaxy and 
a visible quasar in a galaxy w here most star formation has ceased 
(see also lSanders et alj ( Il988h ). At low redshifts, numerous star for- 
mation indicators in quasar host galaxies are seen, e.g. high IR lu- 
minosities, PAH emission, strong (sub)-mm emission and strong 
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CO em is sion (iHao et alj|2005l : ISchweitzer et alj|2006l : iFarrah et al] 
l2007al lbl: iNetzer et al.ll2007ir At high redshifts (z ~ 2-6), obser- 
vations at sub-mm and mm wavelengths (rest-frame far-infrared 
(FIR) to sub-mm, depending on wavelength and redshift range) 
of optically luminous quasars suggest that up to 30% of quasars 
als o fall with in the ULIRG ran ge, similar to the IR quasars s tudied 
bvlCaoet alJfc oOS) ( Carilli e t ai]|200ll: lOmont et al.ll200ll, l2003l: 
ICox et al.ll200a:lBeelen et al.ll2006l : lHao et alfeOOShjcoppin et all 
( l2008h compare dynamical, gas and SMBH masses of ten z ~ 
2 sub-mm detected quasars and z ~ 2 sub-mm galaxies (SMGs), 
which are high redshift counterparts to ULIRGs, though less ex- 
treme with more evenly distributed star formation instead of a lo- 
calised starburst jMenendez-Delmestre et al ]|2009'). They find that 
the fainter half of their quasar sample could likely be 'transition ob- 
jects' between SMGs and luminous quasars, based on their SMG- 
like surface densities and their proximity to the local Mbh/MspJi 
relation, since luminous quasars tend to lie above this relation and 
typical SMGs tend to lie below. The highest redshift sub-mm- 
selected source currently known has also been observed to posses 
similar stellar and gas masses to this z ~ 2 sample of transitio n 
objects, as well as a small AGN contribution jCoDpin et al1l2009l) . 
These quasars eventually may produce enough energy from accre- 
tion to effectively blow out the surrounding gas, halting star for- 
mation and clearing the vi ew to an optically bright quasar before 
halting their own growth jWvithe & Loebl I2OO3I : lOi Matteo et alj 
I2OO5I : iHopkins et al.ll2008l) . This regulation of star formation by 
the accreting SMBH could naturally produce the observed black 
hole-galaxy mass correlation, and different stages of the process 
should present different SFRs, with initially high SFRs declining 
as the AGN becomes more dominant, culminating in a phase of 
strong SMBH accretion and little or no ongoing star formation 
( iKauffmann & Haehnellll2000l : lGranato et al.ll2004l) . 

FIR luminosities can be used to determine the star forma- 
tion rates in quasar host galaxies. The FIR emission is due to 
dust heated eith er by star for mation or q uasar emission (see 
iHaas et alj j2003l) for discussion). [Lutz et alj j2007ll2008l) find ev- 
idence that FIR emission in quasar host galaxies at all redshifts 
is caused by dust heated by star formation and not the quasar. 
The strength of PAH emission, which is found almost exclusively 
in star forming regions, tightly corresponds to the strength of 
FIR emission in the s a me hosts 1 Dale et all I2OOII : ICalzetti et al.l 
I2OO7I : I Lutz et 'ZI I2OO8I) . I Beelen et al. measure FIR emis- 

sion from six high redshift quasars and derive FIR to radio spec- 
tral indexes consistent with local star forming galaxies without 
AGN. The FIR emission is seemingly uncontaminated by hot- 
ter dust potentially heated by the AGN. These results along wit h 
others (see for example, Efstathiou & Rowan-Robinsonl ( Il995l) : 
ISerieant & Hatziminaogloj ( l2009h ) further substantiate the claim 
that the FIR luminosity in quasars is dominated by star formation 
and not by the AGN. 

FIR luminosity traces ongoing star formation, but past star 
formation also can be observed indirectly by measuring chemical 
abundances. Several studies have found that high-redshift quasars 
typically have metallicities greater than or equal to solar metallic- 
ity in the broad emission line region (BL R), which re quires signifi- 
cant previous star formation in the host dHamann & F erland 1999; 
iDietrich et al.ll2003l : IWamer et al.ll2004 fNagao et alj|2006) . These 
BLR studies rely on emission line ratios such as N V/C IV and 
Si IV+O IV]/C IV that are sensitive to the selective enrichment of 
nitrogen as a secondary element and/or to the decreasing tempera- 
tures and increasing metal-line saturations that occur in metal-rich 
gas. The main result for typically solar or higher metallicities near 



quasars has bee n corroborat ed by independent studies of the narrow 
emission lines ^Groves eta l. 2006; Nagao et al. 2006) and narrow 
absorption lines (Hamann & Ferland 1999; D'Odorico et al.ll2004l ; 
iGabel et af> 2006; Simon & Hamann 201CD in quasar spectra. The 
metal-rich BLR result is true even for the highest redshifts stud- 



ied, e .g. iPentericci eTZI ( l2002h ; Ijiang et al.l j2007h ; ljuarez et al] 
( l2009l) . with redshifts out to z = 6.4. There is no kno wn change 
in metallicity with redshift (see iMatsuoka et"ai] ( 12009' ) and refer- 
ences above). Simple chemical evolution models for quasars and 
elliptical galaxies find that galactic centres tend to be more metal- 
rich than their halos, and a centralised starburst can enrich the 
galactic centre to superso l ar abundances in a s hort ti me (^10^ yr) 



(Friaca & Terlevich 1998; Hamann & Ferland 1999; Grana to et al 
| 2001 ; Hamann et al. 2002 ; Granato et al, .2004 ; Hamann et al 
booi lJuarez et al.ll2009l) . These models, combined with the con- 
sistently super-solar gas abundances observed in quasar environ- 
ments, imply that quasars tend to emerge after or near the end of 
(potentially) short centralised star formation epochs. If the quasar 
phase emerges when star formation is on the decline, less ad- 
vanced environments might ha ve higher SFRs and lower abun- 
dances (jGeorgakakis et al. 



ght ha y 

Ullooi) 



The absorption lines in quasar spectra provide additional in- 
formation about outflows that might be related to the blowout of 
gas from the host galaxies, and perhaps, about the gaseous rem- 
nants of recent galaxy mergers. These absorption lines may be 
more common in quasars with higher SFRs if merg ers and inter- 
actions tri gger the s t ar formation as in ULI R Gs iWevmann et al.l 
1991; Be cker et al.l I2OO0I ; iRichardj I2OOII ; iRupke et alj |2005|; 
(Georgakakis et alj 1200 91)" Quasar 'associated' C IV absorption 
lines (AALs), near the emission redshifts with velocity widths 
less than 500 km s^'and broad absorption lines (BALs), with ve- 
locity widths greater than a few thousand kms^', are examples 
of these potential merger, interaction and outflow signatures. In 
systems with recent mergers, remnants from the interaction may 
manifest as a greater incidence of low-velocity (< 2000 km s^') 
AALs. About 25% of bright quasars contain AALs with rest equiv- 
alent width s (REWs) of >0.3 A, and around 10% of quasars 
have BALs jGangulv etaj]|200ll;IVestergaarjl2003l;lMisawa et al.l 



2003*; 'Trump et al.' '2006*; 'Nestor et al.' '2008*; jwild et al.l |2008|; 
Gibson et al. 2008; Rodriguez Hidalgo et al. 20l3). A higher in- 
cidence of AALs or BALs may occur in the quasars with higher 
SFRs if these quasars have more recently experienced an interac- 
tion and/or there is a progression in quasar outflow characteristics 
with time. 

In this paper, we present an exploratory observational study 
examining whether SFR in the host galaxies correlates with metal- 
licity in the near-quasar environment. We measure metallicity in 
the quasar BLR from the rest frame ultraviolet (UV) spectrum and 
estimate galactic SFRs from the FIR luminosities. 

The data and analysis are described in § |2] and § [3] The re- 
sults and discussion are presented in § |4l with a brief summary in 
§[5] We adopt cosmological parameters Hq = 70 km s^' Mpc^', 
Q.m = 0.3, and Q.\ = 0.7 throughout this work. 



2 DATA 

We select quasars for this study from t he sample o b served 
wit h MAMBO a t IRAM at 1.2 mm by ICarilli et alj ( 1200 ll) 
and 'Omont et al.' ("2001', '2003") and SCUBA a t JC MT at 850nm 
by Isaaketal. (2002), McMah on et alj d 19991) and IPriddev et all 
( I2003i) and all compiled bv iHao et al. ( 2008[) . The quasars all were 
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selected to be optically bright with absolute B-band Magnitude 
Mb < -26.1 for the CariUi et al. objects, Mb < -27.0 for the 
Omori t et al. object s and M b < -27. for t he 'McMa hon et al.1 
(Il999h: llsaak et al.l ( l2002h and IPriddev et al . (2003) objects. 
ICarillietal.U200lh observed a representative sample of 41 out of 
the more than 100 quasars with z ^ 3.6 found as part of the Sloan 
Digital Sky Su rvey (SPSS ) Galactic Ca p and Southern Equatorial 
Stripe survey. Omo ntetal.l l l200lLl2003h observed a random selec- 
tion of 97 radio quiet, optically luminous sources from the multi- 
colour Palomar Digital Sky Survey a vailable from G. Djorgovski's 
web pag43with 3.9 < z < 4.5 a nd the ' Veron-C etty & Veronl fioooh 
catalogue with 1.8 < z < 2.8. iMcM ahon et al.l il99'9r selected a 
small sample of 6 bright z > 4 radio quiet quasars from the APM 
survey tStorrie-Lombardi et al .Il996l) . llsaaketalj l l2002l) selected a 
larger sample of the 76 most UV lumi nous z ^ 4 rad i o quie t quasars 
known at the time of observation, and lPriddev et al.l j2003h selected 
a complimentary sample of 57 z ^ 2 quasars from various large sur- 
veys. 

We cross-reference this sample with the optical quasar spec- 
tra in the SDSS data-release 6, finding 116 objects with available 
spectra. The SDSS spectra from data release 6 have resolution 
R= X/5?i ~2000 and wavelength coverage 1 = 3800 - 9200 A 
jAdelman-McCarthv & et al.l [20081) . We consider only those ob- 
jects with redshifts between 2.17 and 4.75, compatible with SDSS 
spectral coverage of the full 1200 - 1600 A rest-frame wave- 
length range. Spectra missing regions within the specified wave- 
length range are further excluded from the analysis. Five spectra 
with very low signal to noise ratios (S/N) also are rejected. The 
final sample consists of 34 optical SDSS spectra with a range of 
sub-mm brightnesses and a redshift range of 2.2 < z < 4.6. 

The MAMBO 1.2 mm and SCUBA 850 ^m observations 
correspond to 160 - 400 /Jm rest wavelengt hs, depend i ng on 
quasar redshift. We follow the prescription in iHao et alj feOOSh 
to convert from observed flux into FIR luminosity at 60 ^m, 
Lgo = XLx{60 flm), in which we assume that a grey-body 
spectrum describes the rest frame FIR spectral energy distribution 
(SED) with a dust te mperature of 41 K and dust emissivity index 
of 1.95, as detailed inlPriddev & McMahonI (1200 ih . We find agree- 
ment with lHao etal] ( l2008 l) to within 10%. 

The Lgo, absolute B-band magnitude (M^) and bolometric lu- 
minosity (Lboi) for each object are shown in Figure[T] and listed in 
Table[T] The FIR luminosity, log(L6o/L0), ranges from ^12.1 to 
13.4, is dominated by star formation and falls roughly within the 
ULIRG range. We estimate the star formation rate of each quasar 
host from L50, corrected to exclude the small quasar contribution 
by assuming the hosts follow the same regress ion line as typica l 
quasars for Lboi vs. Lgo as shown in Figure 1 of iHao et alj j2008l) . 
This corrected Lgo is then used in Hao et al.'s equatio n 2, which 
is der ived from the Kennicutt star formation rate law (Kennicutt 
Il998l) . and the estimated star formation rates for each object are 
also listed in Table [T] . We convert the Mb from the observation 
papers to our adopted cosmology, and obtain spanning the 
range of -26.7 to -29.4, which corresponds to log(L|,oi) of 47.2 
to 48.3 log(erg s^'). Lboi measures the quasar black hole accre- 
tion luminosity, and is estimated using a bolometric correction fac- 
tor of 9.74 applied to the mon ochromatic continu um luminosity 
AL;t(4400 A) (~ Mb) following |\^stergaardl ( l2004h . 



' They use an Einstein deSitter cosmology, Q.m = L^a = OMo = 
50 kins''. 

^ http://astro.caltech.edurgeorge/z4.qsos 
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Figure 1. Quasar FIR luminosities, absolute B magnitude and bolomet- 
ric luminosities. Horizontal dashed lines represent the division between the 
three FIR luminosity bins, crosses are FIR-detected quasars and arrows are 
upper limits. Filled triangles are average luminosities for each FIR luminos- 
ity bin. 



We separate the sample into three bins based on Lgo such 
that each bin contains a similar number of objects: FIR bright 
quasars with log(L6o/LQ) ^ 13.17, FIR intermediate quasars 
with 12.8 < logCLeo/L©) ^ 13.1 and FIR faint quasars with 
log(L6()/L0) < 12.75, in which the FIR faint luminosities are upper 
limits. Any objects with logCLgo/L©) upper limits above 12.75 were 
not included in the study. The FIR-bright and FIR-faint bins each 
span the redshift range from ~2.2 to --^4.4. The FIR bright bin also 
includes one object with z --^4.6, and the FIR intermediate bin spans 
a smaller redshift range from 3.7 to 4.4. The divisions between the 
Lgo bins are represented in Figure[T]by horizontal dashed lines and 
average values in each bin for Lboi and Lgo are denoted by filled 
triangles. These average values are listed in Table|2]for each lumi- 
nosity bin in columns 2 and 3, along with the corresponding star 
formation rates in column 4, and several abundance indicators, dis- 
cussed in ij's |3.3| and|4l 



3 ANALYSIS 

3.1 Composite Spectra 

To make comparisons between different L^y, we create one com- 
posite spectrum from the SDSS spectra comprising each of the 
three FIR bins in the sample described in § |2] To create the com- 
posites, we begin by shifting each spe ctrum to its rest wavele ngth 
using the redshifts provided by SDSS dSubbaRao eTaLll2002h . We 
manually inspect each individual spectrum and remove absorption 
features by interpolating across the affected regions. The presence 
of narrow absorption lines and their relationship to Lgo is discussed 
in i) |3.4| below. Then we average together the spectra in each bin to 
create the final composites. To ensure that no single spectrum is 
dominating the composite we also compute the median spectra for 
each bin, compare them to the average spectra used throughout the 
rest of the analysis, and confirm that the two composite types are 
well-matched. 

Our analysis is limited to the spectral region between 
Lya I2I6 A and C IV 1550 A. These limits are imposed by the in- 
complete spectral coverage at longer wavelengths and by suppres- 
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Table 1. Quasar Sample 





z 


logLboi 


Ms 


logLeOfim 


SFR 














M0yr 


Ref.° 


FIR Bright 














PSS J1057+4555 


4.12 


48.1 


-29.0 


13.2 


3720 


OOl 


PSS J1347+4956 


4.56 


47.9 


-28.5 


13.3 


4690 


OOl 


J144758.46-0()5055.4 


3.80 


47.2 


-26.8 


13.3 


6350 


COl 


J154359.3+535903 


2.37 


47.9 


-28.4 


13.3 


5550 


O03 


PSS J1248+3110 


4.35 


47.6 


-27.8 


13.3 


6220 


OOl 


PSS J1418+4449 


4.28 


48.0 


-28.8 


13.3 


5320 


OOl 


PSS J0808+5215 


4.45 


48.1 


-28.9 


13.3 


5350 


OOl 


Jl 10610.8+640008 


2.19 


48.3 


-29.4 


13.4 


4720 


O03 




2.63 


48.0 


-28.6 


13.2 


3290 


P03 


J142647.82+002740.4 


3.69 


47.3 


-26.8 


13.2 


4580 


COl 


FIR Intermediate 














J02323 1.40-0000 10.7 


3.81 


47.2 


-26.7 


12.8 


1870 


COl 


PSS J1535+2943 


3.99 


47.4 


-27.3 


12.8 


1710 


OOl 


J012403.78+004432.7 


3.81 


47.9 


-28.4 


12.9 


1230 


COl 


J0255 18.58+004847.6 


3.97 


47.6 


-27.8 


12.9 


1670 


COl 


J015048.83+004126.2 


3.67 


47.7 


-27.9 


12.9 


1920 


COl 


J141332.35-004909.7 


4.14 


47.4 


-27.3 


12.9 


2330 


COl 


J025 112.44-005208.2 


3.78 


47.3 


-26.9 


13.0 


2580 


COl 


Jl 11246.30+004957.5 


3.92 


47.6 


-27.7 


13.0 


2530 


COl 


PSS J1317+3531 


4.36 


47.6 


-27.7 


13.1 


3330 


OOl 


FTP Fnint 














TOI^SIQ 61-000910 1 


4.20 


47.5 


-27.4 


5^12.5 


^430 


COl 


T140SS4 07-0000^7 


3 55 


47 5 


-27.6 


^12 7 


^^850 


COl 


Mi^AA Q 9'^ ooni n 
J zz j'+ 1 y. z j-uuui J J. u 


j.Oo 




lf\ Q 

-zo.y 


^19 1 

5^ IZ. 1 






PSS J1443+5856 


4.27 


48.0 


-28.7 


^12.4 


~0 


OOl 


J225759.67+001645.7 


3.75 


47.5 


-27.4 


^12.2 


-.0 


COl 


HS B0808+1218 


2.26 


47.6 


-27.8 


^12.4 


-.0 


P03 


HS Bill 1+4033 


2.18 


47.7 


-27.9 


^12.7 


^550 


P03 


LBQSB1210+1731 


2.54 


47.8 


-28.1 


^12.4 


~0 


P03 


LBQSB1334-0033 


2.80 


47.6 


-27.8 


^12.6 


^610 


P03 


BR 1600+0729 


4.35 


47.8 


-28.3 


^12.7 


^280 


M99 


PSS J0852+5045 


4.2 


47.7 


-28.0 


^12.4 


~0 


102 


PSS J0957+3308 


4.25 


47.8 


-28.1 


^12.7 


^650 


102 


PSSJ1618+4125 


4.21 


47.4 


-27.2 


^12.7 


^990 


OOl 


PSS J1633+1411 


4.35 


47.6 


-27.6 


^12.6 


^570 


OOl 


J142329.98+004138.4 


3.76 


47.2 


-26.6 


^12.2 


^180 


COl 



"OOl lOmont et alj i200ll). O0 3 lOmont et alj j2003h 102 llsaak et alj \200% . COl ICarilli et alj 
j200lh . P03 |Priddev et alj J2003h . M99 lMcMahon et alj i\9^ 



Table 2. Metallicity from Emission Line Flux Ratios 



Composite spectra 


<log Lboi > 
log(erg s"') 


<log Leo > 
log(erg s"') 


<SFR> 
M^yr-' 


Flux Ratio 
N V/C IV 


Z/Z0 
N V/C IV 


Flux Ratio 
Si IV+O IV]/C IV 


Si IV+O IV]/C IV 


FIR Faint 


47.60 


^46.07 


^340 


1.2 


7.9 


0.3 


3.0 


FIR Int. 


47.52 


46.52 


2130 


1.4 


9.6 


0.3 


3.0 


FIR Bright 


47.89 


46.87 


4980 


1.5 


10.9 


0.4 


6.5 



sion in the Lya forest starting in the blue wing of the Lya emis- 
sion line and extending to shorter wavelengths. Before normahsing 
the composite spectra, we visually compare the continuum slopes 
of the three composites. We find no significant trend for redden- 
ing with FIR luminosity in the continuum slopes of the composite 
spectra. Differences between the composite slopes are negligible 



compared to the dispersion of slopes among the individual quasars 
making up each composite. 

We fit a power-law continuum to the FIR-intermediate com- 
posite spectrum using wavelength regions devoid of absorption or 
emission, 25 A wid e and centred at A 1460, and A 1770, following 
IWamer et all j2003l) . Unfortunately, the FIR-faint and FIR-bright 
composites have incomplete spectral coverage past ~A1700 A, so 
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we fit the continuum using tlie available wavelength region around 
X 1 460 A, plus t he reg ion from A 1335-1355 A, which is also used 
bv I Juarez et alJ ( 1200 9') . and which we determine by visual inspec- 
tion to be devoid of absorption or any detectable emission. The 
normalised average composite spectra are shown in Figure[2] 



3.2 Emission Line Flux Ratios 

We measure several quasar broad emission line fluxes to cal- 
culate line flux ratios, which we then use to estimate the 
gas phase metallicity in the near-quasar environment. Flux 
ratios are meas ured for the emissi on lines that have been 
shown by e.g. iHamannetal] ( l2002h : ljuarez et al] j2009l) . to 
give good abundance estimates within the limited wavelength 
range from A 1216 A to A 1600 A: N V 1 1240/C IV A 1549 
and Si IV Al397-l-OlV] Al402/ClV A 1549. These emis- 
sion lines are labelled in Figure |2l Other ratios such as 
N III] A1750/O III] A 1664, Nv Al240/He II A 1640 or 
N III A991/C III A977 are too weak given the S/N in the spectrum 
and/or they fall in a poorly characterised region of the spectrum. 
The measured line flux ratios are listed in columns 5 and 7 of Ta- 
ble[2l 

Because the N V emission is strongly blended with the 
Lya emission, which is itself significantly degraded by absorption 
in the Lya forest, C IV is the only strong, non-blended emission 
line in the spectra. Thus, the C IV emission lines are fit with gaus- 
sians and used as templates for the other emission lines. The fits to 
C IV and the scaled fits to Lya and N V are shown for the three 
composites in Figure [3] The fits are performed using GATOR- 
PLOT, an IDL program written by C. WameJ^, and are carried 
out using the smallest number of gaussians possible to provide a 
good match to the data: one broad and one narrow component for 
each of the two doublet lines at A 1548 and A 1551 A. The two gaus- 
sians comprising the A 1548 emission line fit are held at fixed full 
width at half maximum (FWHM) and are scaled in total flux at the 
fixed central wavelength of the Lya A 1216 A and N V A 1239, 
1243 A doublet emission to match these line strengths. Because 
C IV is known to be blue-shifted relative to Lya an d the nominal 
quasar redshift by ~ 3 10 km s^' jTvtler & Fanll993) . the locations 
of the N V and Lya centroids are shifted from their laboratory val- 
ues by this amount relative to C IV. 

To accurately determine the N V flux, the red Lya wing be- 
neath the N V emission must be characterised. To do this precisely, 
we scale the C IV fit to align with the non-absorbed region of Lya, 
which in some cases is quite small 20 A). The scaled fit may 
even rise above the data where heavy Lya forest absorption has 
eaten into the emission, as is the case for all of the three compos- 
ites, or miss the peak of Lya as is the case in the FIR faint com- 
posite; both cases are clearly shown in Figure[3] In order to better 
match the Lya emission peak in the FIR faint composite, we allow 
the Lya gaussian to be ~ 15% narrower than the C IV gaussian. 
Narrowing the gaussian for the Lya emission line does not have a 
noticeable effect on the resulting N V strength. Regardless, a pre- 
cise match to the Lya emission line is not required for this anal- 
ysis, and does not effect the final N V emission results (see e.g. 
[Baldwin et al.l (2003) for further discussion). 

The C IV and N V doublets are fit at their respective fixed 
separations and with a 3:2 intensity ratio (halfway between the al- 
lowed 2: 1 and 1 : 1 intensity ratios) for the shorter and longer wave- 



length respectively, as in [Baldwin et id] ( l2003h . Possible broaden- 
ing in N V is predicted if the N V emission for ms at a larger 
distance from the quasar than the C IV emission. [Peterson et al.l 
(2004) and Peterson (2008) find that N V forms twice as close to 
the quasar as C IV, so if the line widths are controlled by virial 
motions, N V could be up to ~ \/2 times broader than C IV. 
Lya, on the other hand, forms at about the same distance from 
the quasar as C IV and should be no broader (an d probably nar- 
rower) than C IV ( |Petersonetlj]|2004l : [Peterson 2008). There- 
fore, the N V profiles we adopt, which have the same FWHM's 
as C IV, may be narrower than the actual N V FWHM's, while the 
Lya FWHM could be slightly broader. We estimate the uncertain- 
ties in the N V line strengths due to uncertainty in the continuum 
placement by repeating the line fits with different reasonable con- 
tinuum heights. We also include estimates of N V line strength 
uncertainties due to blending with Lya and the poorly defined 
N V FWHM by repeating the line fits with a range of different 
N V and Lya FWHM's. The overall uncertainties, combining 
these effects, are less than a factor of 1.5, so that the measurement 
errors in the N V/C IV ratios are also less than a factor of 1.5. 

The Si IV-l-0 IV] emission line is weaker than the C IV and 
N V emission lines in the composite spectra, however it also 
is isolated from other emission. This isolation means that mea- 
suring the equivalent width of the emission without fitting, us- 
ing IRAlj^ produces accurate fluxes. We confirm the accuracy of 
the equivalent width measurements by comparing the previously 
measured fluxes from gaussian fits to C IV to the fluxes from 
C IV equivalent widths, and find compatible results. We mea- 
sure Si IV-l-O IV] and C IV equivalent widths to calculate the 
Si IV-l-O IV]/C IV line flux ratios listed in Table|2] To estimate un- 
certainties in the Si IV-l-O IV] line strengths, we vary the equivalent 
width measurements by using a range of continuum heights and 
wavelength cutoffs for the edge of the emission lines. The overall 
uncertainties in the Si IV-l-O IV] line strengths are less than a factor 
of 1.2, so that the measurement errors in the Si IV-l-O IV]/C IV ra- 
tios are also less than a factor of 1.2. 



3.3 Metallicity 

We convert the N V/C IV flux ratios into metallicities using 
the theoretical relationship in which nitrogen abundance in- 
creases relative to carbon as metallicity increases because of sec- 
ondary enrichment (CNO nucleosynthesis) processes, and deter- 
mine the average metallic i ty for each Lgo composite ( Shields 197^; 
Hamann & Ferland 1 9931 , [I999I) . The correlation is characterised 
by Ha mann et aL C2002[) with recent updates to the solar abundance 
ratios by Dhanda et alj ( [2007h . We find supersolar metallicities for 
all three composite spectra, as listed in column 6 of Table|2] 

Similar characterisations for the Si IV-l-O IV]/C IV flux ratio- 



Nagao et al. Ll2006^, using 
Hamann et al.) C2002) . We 



metallicity correlation are performed by 
the older solar abundances also used by 
apply the latest corrections for the solar abundance ratios for this 
characterisation and determine that the metallicities are also super- 
solar according to this ratio, as shown in column 8 of Table|2] 

There is no significant trend among the three composites in 
the line ratios, based both on our measurements of these ratios and 
on a visual inspection of the stacked composites in Figure (2] and 



^ http://www.astro.ufl.edur wamer/GatorPlot/ 



^ IRAF is the Image Reduction and Analysis Facility, supported by NOAO 
and AURA Inc. 
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Figure 2. Nonnalised composite spectra for FIR-faint (solid line), FIR-intermediate (dashed line) and FIR-bright (dotted line) for the (rest) wavelength range 
from Lya to C IV. Prominent emission features are labelled. The continuum has been normalised to 1, as described in the text. 



correspondingly there is no significant trend in metallicity. The av- 
erage metallicities across the three composites are Z'-^9.5 Zp, for 
the N V/C IV ratio and 4.2 Z© for the Si IV+0 IV]/C IV ratio. 



We note that the metalMcities inferred from the two line ratios 
can differ by as much as a factor of '-^2 in the same spectrum. This is 
typical of the dispersion found between different hne ratios in other 
studies, and it might be a good indication of the theoretical uncer- 
tainties jPietrich et a l. 2003; Nagao et al. 2006; Hamann & Simon 
I2OO9 !). We compare our metallicities to the metallicities in the 
Lboi range 10"*^ to 10"*** erg s ' sampled by a large emission line 
study bv lWamer et alj ( |2004) . who find metallicities of Z --^3-6 Zq, 
broadly consistent with the metallicities found in this sample. We 
also note that the results based on N V/C IV(and N V/He II) tend 
to be higher than other line ratios, and thus our best guess at the 
metallicity in the current sample overall would be near the lower 
end of the range 4- 9 Z© ( iHamann et alj2002l ; [Baldwin et ahlioosl ; 
iNagao et al.l l2006h . The relative metallicities between the three 
composites, which smooth over object-to-object scatter, are robust 
for differences greater than a factor of 2-3, based on the uncertain- 
ties in the line ratio measurements plus the theoretical uncertainties, 
and therefore useful for spotting strong metallicity trends among 
the different Lso bins. 



3.4 Absorption Lines 

We inspect each SDSS spectrum used to make the composites and 
count the number of C IV AALs with REW > 0.3 A per spec- 
trum. We find that 17/34 quasars have at least one AAL and 6/34 
have absorption too broad to be classified as an AAL, yet too 
narrow to be classified as a BAL, with widths more than twice 
the thermal width and velocity dispersion expected for unrelated 
gas in the line of sight, which is indicative of g as forming in 
an outflow l lBahcall et al]|l967l ; [Young et al.lll982l) . These statis- 



numbers involved in this study 1 


Vestergaard 


I2OO3I; iTrumc et al.l 


20061; Nestor etalj 20081; Wildetal. 20081; 


Gibson etal. 20081; 


Rodriguez Hidalgo et al.l2010'). AALs and outflow lines acoear in 



each of the three Lgo bins, with no significant differences in their 
occurrence fractions among the bins that would indicate a trend 
with Lgo, particularly given the small number of objects in each 
bin. 



4 DISCUSSION 

The average line flux ratios N V/C IV ~ L35 and 
Si IV-l-O IV]/C IV ~ 0.33 across the three composites, as 
shown in columns 5 and 7 of Table [2] correspond to average 
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metallicities of Z ~ 9.5 and 4.2 Zr;i respectively, shown in columns 
6 and 8 of Table (2] There is no significant trend in metallicity 
with Lgo in this sample. We note that the FIR bright bin has the 
most luminous average L(,oi, while the FIR faint and interme- 
diate bins have roughly equal, and less luminous average L(,oi. 
This difference in l^bo\ could appear in the relative metallicities 
among the bins because more luminou s quasars tend to be more 
metal-rich than less luminous quasars ( Hamann & Ferland 1999'; 
I Warner e~ al. 2004; She mmer et al. 2004). However, the absence of 
a corresponding gradient in the metallicity is not surprising given 
the relatively high average of the total Lboi, the range in average 
Lboi among the composites of less than 0.4 dex (Table which 
corresponds t o a ra nge in metallicity of no more than 0.2 dex in 
IWamer et al. I j2004h . and the sensitivity of the data to trends in 
metallicity greater than --^0.3 dex between the three L5() bins. 

The fact that the data show no significant metallicity evolution 
with changing SFR (Lgo) could be due to the small sample sizes 
being affected by object-to-object scatter or by the small dynamic 
range in L50, with values that are all within the luminosity/SFR 
range of powerful ULIRGs. 

Absorption features such as AALs and BALs and reddened 
continua may be more numerous among quasar environments that 
have undergone a recent merger or are particip ating in the b l owout 
of gas that is revealing the visible quasar teecker et ai] I2OOOI ; 
iRichardsl I2OO1I ; iD'Odorico et alj |2004| ; iGeorgakakis et al.l l2009h . 
Although these are optically selected quasars chosen to be bright 
in the rest-frame UV and strongly reddened sources would be ex- 
cluded, there could still be a more subtle trend in sub-samples by 
LgQ. We do not find a trend in the number of AALs per quasar or in 
the continuum slope per L50 bin. Each bin appears to have roughly 
the same percentage of quasars with AALs. The absence of a trend 
could be due simply to the small numbers in our sample. Or, if the 
consistency across bins is real, it could indicate that there is no pro- 
gression in quasar outflow characteristics with changing SFRs. A 
larger sample size is needed to resolve this question. 

The overall conclusion from this and other studies is that 
quasars appear to be metal-rich at all redshifts, with only 
smal l variations due to SMBH mass and luminosity dependen- 
cies ("Warner et al."2004'; 'Shemmer et ai]|2004l ; iNagao et aij|2006l ; 
Ijiang et a l. 2007; Juarez et al. 2009). The consistently high metal- 
licities suggest that the star formation producing these FIR-bright 
quasars is not the star formation that determines the metallicity of 
the BLR. Instead, the BLR gas must have been enriched prior to 
these FIR-producing star formation episodes, possibly shortly after 
a merger or some other star formation trigger when the starburst 
first began. There are no significant differences among the Lgo bins 
in the outflow fraction or in the amount of debris leftover from 
mergers in the form of absorption, which might correlate with Lgg 
if this is indeed an evolutionary sequence. This scenario is consis- 
tent with models of galaxy and quasar evolution, which predict a 
lack of metallicity evolution. The models also find that quasar ac- 
tivity, though caused by a disruption of some sort, is significantly 
delayed after the disruption event and most of the star formation 
and enrichment in the host is complete even as the quasar phase 
first begins to emerge l Oi Matteo et al.ll2005l ; lHoDkins et al ] |2008l ; 
iLi et aklboOSl) . The FIR-bright quasars may be in some sense less 
mature than the FIR-faint quasars, however, this is not manifest in 
significant metallicity evolution over the timescales probed by this 
sample, which are shorter than a quasar lifetime. 

The reality of quasar-host galaxy formation could also be 
more complicated than any simple monotonic sequence in which 
high SFRs correlate cleanly with lower metallicities and an earlier 



stage of evolution. IVeilleuxetai] ( l2009h suggest a 'softer' ULIRG- 
quasar evolution paradigm with more scatter in the path of an in- 
dividual galaxy's evolution from ULIRG to quasar. They measure 
neon abundances of Z~2.9 Z0 in the nuclear regions of ULIRGs. 
If these high metallicities are ubiquitous in ULIRGs, we will never 
observe metal-poor quasars, regardless of their stage of evolution. 
Other plausible scenarios suggest that AGN and star formation ac - 
tivity in galaxies may be episodic in nature (Di Matteo et al.'2008j). 
Thus, as a quasar emerges at the tail end of a star-forming phase, 
there is no guarantee that its BLR metallicity is linked to this re- 
cent and/or ongoing star formation phase, because several previous 
episodes of star formation and AGN activity could already have 
occurred and enriched the BLR gas to supersolar values. 

An alternative interpret ation is provided by t he wo rks of 
iDave & QppenheimeJ ( |2007|) and iFinlator & Davi l l2008h . They 
postulate a scenario where a star-forming galaxy could quickly 
reach an equilibrium metallicity early in its formation history. 
Higher mass galaxies would have higher abundances, but for a 
given galaxy, the abundance would not change significantly after an 
initial balance of in-flowing metal-poor and out-flowing metal-rich 
gas was reached. The lack of BLR metallicity evolution with star 
formation rate found in this study is consistent with these objects 
having previously reached an equilibrium metallicity, unaffected by 
the observed ongoing star formation. 



5 SUMMARY 

We create composite rest-frame UV spectra for a sample of 
34 quasars observed in the FIR and measure N V/C IV and 
Si IV+O IV]/C IV flux ratios. We convert these flux ratios into a 
metallicity for each of three Lgg bins, which we interpret as SFR 
bins. We find that all three bins have supersolar metallicities of the 
order 4 to 9 times solar, and find no metallicity trend with Lgy. We 
investigate the amount of outflowing gas, as manifest by both AALs 
and somewhat broader absorption lines, in each Lgo bin and find no 
difference in number between the Lgo bins, neither do we find a 
change in continuum slope with changing Lgo, as would be consis- 
tent with more dusty star formation in hosts with higher SFRs. To- 
gether, these results suggest that the ongoing star formation in the 
host is not responsible for the metal enrichment of the BLR gas. 
Instead, the BLR gas must have been enriched before the visible 
quasar phase. We note that we require a substantially larger sample 
size with a broader range in Lgo to confirm that these results are not 
affected by unusual individual objects or biased by the narrow Lgo 
range. 
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Figure 3. Gaussian fits for Lya, and N V in the left panels and C IV in the 
right panels for each normalised composite spectrum. The vertical scale is 
normalised flux. The top panels shows the FIR faint spectrum, the middle 
panels shows the FIR intermediate spectrum and the FIR bright spectrum is 
shown in the bottom panels. The soUd curves are the data and total gaussian 
fits, the dashed curves show the individual gaussian fits for each emission 
Une. The N V and C IV features each have two gaussians, one for each 
doublet member, which are each composed of a narrow and a wide gaussian 
(not shovm). The error spectra for each region are shown as dot-dashed 
Unes. 



